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We take a first step towards the solution of QCD in 1 + 1 dimensions at nonzero density. We 
regularize the theory in the UV by using a lattice and in the IR by putting the theory in a box 
of spatial size L. After fixing to axial gauge we use the coherent states approach to obtain the 
large- A r classical Hamiltonian Ti that describes color neutral quark-antiquark pairs interacting with 
spatial Polyakov loops in the background of baryons. Minimizing TL we get a regularized form of 
the t Hooft equation that depends on the expectation values of the Polyakov loops. Analyzing the 
Independence of this equation we show how volume independence, a la Eguchi and Kawai, emerges 
in the large-iV limit, and how it depends on the expectation values of the Polyakov loops. We 
describe how this independence relies on the realization of translation symmetry, in particular when 
the ground state contains a baryon crystal. Finally, we remark on the implications of our results 
on studying baryon density in large- iV QCD within single-site lattice theories, and on some general 
lessons concerning the way four-dimensional large- N QCD behaves in the presence of baryons. 

PACS numbers: 11.15.Ha,11.15.Pg 



I. INTRODUCTION 



QCD simplifies in the 't Hooft limit of a large number of colors, and as a result it has been a long-standing goal to 
understand the properties of the theory in that limit [l[ , including on the lattice • One alternative to conventional 
large volume simulations is the use the large-TV equivalence of QCD at large volume to QCD with zero volume 
[E EL IE Hi 0, H, H 03 EH (see also the related Ref. [Ilj]). These large- N volume reductions allows one, in principle, 
to study very large values of N ~ 0(100 — 400) with modest resources. 

Volume reduction holds only if the ground states of the large and zero volume theories respect certain symmetries [1[ . 
Unfortunately, in the most interesting case of QCD in four dimensions these symmetries spontaneously break in the 
continuum limit when a naive reduction prescription is used 0, Q • An extension of that prescription is thus required 
and for a recent summary of the literature on this topic we refer to the reviews in the recent Refs. (l3l ITU]. 

In the case of two space-time dimensions - the 't Hooft model - a naive large- TV volume reduction is expected 
to hold and so this theory is generally thought to be completely independent of its volume. In the current paper 
we analyze this volume dependence. Our motivation is two-fold. Firstly, the 't Hooft model is analytically soluble 
at large-TV. Thus we can explicitly see how large-iV volume reduction works in this case, and what may cause it 
to fail. This topic was also addressed for zero baryon number by the authors of Ref. [HI, and our treatment here 
differs from that paper by being manifestly gauge invariant, by going beyond zero baryon number, and by using the 
lattice regularization. Our approach also makes a direct connection with Eguchi-Kawai reduction, and shows how the 
expectation values of spatial Polyakov loops play a crucial role in the validity of volume independence. 

Secondly, this paper is a prelude to our companion publication Ref. [l|| where we use the formalism presented here 
to solve the theory in the presence of nonzero baryon density. Considering the current incomplete understanding of the 
way four-dimensional QCD behaves at low temperatures and large (but not asymptotic) baryon densities, we believe 
that such a study is useful. Also, there exist certain confusions in the literature about the way large- N gauge theories 
behave at nonzero baryon number [13], and seeing how these confusions go away in the soluble two-dimensional case 
is very helpful. 

Surprisingly, QCD in two dimensions and nonzero density has not been solved yet : While Ref. [l8[ studied only one 
and two baryons in an infinite volume, then Ref. [li^ attempted to extend this but restricted to either (1) translational 
invariant states which were seen to be inconsistent, or (2) a particular translational non- invariant ansatz for a baryon 
crystal in the vicinity of the chiral limit. Since 1 + 1-dimensional baryons become massless for massless quarks, it is 
natural to expect that they behave very differently than four-dimensional massive baryons. Furthermore, most of the 
current literature on the 't Hooft model has so far focused on its infinite volume limit where a certain set of gluonic 
zero modes are irrelevant. With a finite density of baryons, however, these become important and cannot be neglected 
(at least if the density is increased by fixing the baryon number and decreasing the volume). Thus, given the current 
status surveyed above, it seems wise to study the dense 't Hooft model for arbitrary quark mass, by making as few 
assumptions on the form of the ground state as possible, and by incorporating correctly the gluonic zero modes. In 
this paper we develop the machinery to achieve this goal. For the actual solution of the theory for arbitrary baryon 
numbers we refer to Ref. [l6j and for all other discussions on the way nonzero chemical potential affects the system 
to Ref. [II. 

Former studies of the 't Hooft model used a plethora of mathematical methods - for example see [2l[ for some 
papers relevant to this work. Common to all these is the need to control the severe IR divergences of this two- 
dimensional model. A particular clear approach, that we will follow in our study, is the one advocated in the seminal 
Refs. [22|, There, one works in the Hamiltonian formalism defined in a spatial box of side L, and uses the axial 
gauge to remove all redundant degrees of freedom. This approach is also most suitable for our purpose of investigating 
the L-dependence of this Hamiltonian's ground state. 

The outline of the paper is as follows. In Section [II] we present the details of the Hamiltonian approach to lattice 
QCD. A reader who is familiar with this approach can skip to Section IIIII where, by generalizing Refs. [22l . [23j to 
the lattice, we show how to fix the axial gauge in the Hamiltonian formalism. Since such gauge fixing is less familiar 
than the gauge fixing in the Euclidean formalism, we do so in detail. Next, in Section llVl we show how to resolve 
Gauss law and re-write the electric fields in terms of the fermion color charge densities. This rewriting can be done 
for all components of the electric field except for those conjugate to the eigenvalues of the spatial Polyakov loops. 
This set of eigenvalues and their conjugate electric fields is what we refer to above as zero modes, and in Section [V] 
we focus on them. Specifically, we show how to represent the zero modes' electric fields in the Schroedinger picture 
as differential operators. The end result of Sections HT1IVII is a Hamiltonian that depends only on the fermions and 
on the zero-modes, with an overall color neutrality enforced on its Hilbert space. For the convenience of the reader 
we summarize this emerging structure in Section IVII We then turn to find the ground state of this Hamiltonian. At 
large- TV this is done in two steps : (1) Solution of the gluon zero modes dynamics - discussed in Section fVII Al (2) 
Solution of the fermion sector - Section IVII Bl In the second step we use the coherent state approach of Refs. [25| 
which seems particularly suitable for our problem. The end product is a regularized form of the 't Hooft classical 
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Hamiltonian describing color neutral operators that correspond to quark-antiquark pairs and Polyakov-loops (that 
wrap around the spatial circle), and that interact in the presence of a fixed overall baryon number. 1 In Section [Villi 
we survey other relevant works in the literature that obtain a similar Hamiltonian, pointing out the way they differ 
from our approach. In Section IIXI we analyze the resulting Hamiltonian and its L dependence for arbitrary baryon 
number B. We show how large- N volume dependence emerges and that for it to hold we need to assume that the 
ground state has some degree of translation invariance. We also show how it can be violated by giving the Polyakov 
loops nonzero expectation values. 2 An interesting phenomena occurs when the ground state contains a baryon crystal 
and we show how a 'soft' form of volume independence emerges. This leads us to remark on the way our results affect 
studies of nonzero chemical potential that try to rely on large- iV volume independence. We conclude in Section IXl by 
noting some general lessons one can learn about the way large- N QCD behaves in the presence of baryons. 



II. HAMILTONIAN QCD IN 1 + 1 DIMENSIONS : A BRIEF REMINDER 

In this section we introduce the Hamiltonian formalism of lattice QCD restricted to one spatial dimension and one 
flavor. A reader familiar with this formalism can skip to Section [IIII 

This Hamiltonian of lattice QCD was first introduced by Kogut and Susskind in 1975 [26| . shortly after Wilson's 
Euclidean formulation. This canonical formalism defines the theory of strong interactions on a spatial lattice with 
lattice spacing a and continuous time t. In one dimension a lattice site is denoted by a single index x taking integer 
values. We also use x to denote the lattice link that is to the right of the site x. Because we define the theory on a 
finite box we set x — 1,2,..., L s , where L s = L/a is the number of lattice sites, and L the physical length of the box. 
The boundary conditions of the gauge fields are taken to be periodic. In this section and throughout the paper we 
use standard lattice notations and so the factors of the lattice spacings are implicit and all fields are dimensionless. 

The quantum fields that describe quarks are the fermion fields ip% a that live on the sites x of the lattice. They have 
color indices a = 1, . . . , N, and Dirac indices a = 1, 2 (again, recall we are in 1 + 1 dimensions). The Fermi fields 
obey the following anti-commutation relations 

{vr,^ b/3 } = <www- (i) 

Choosing to work in the temporal gauge that fixes Aq = removes one degree of freedom (and its conjugate 
momentum). This leaves the spatial gauge field operators (U x ) ab living on the lattice link between x and x + 1. For 
definiteness we note that the operator ([/,£) , is defined to be equal to ((U x ) ba ) , and so the following holds as a 
operator identity 

E ((^U 1 (^)ac = E (^L (^)ac = S bc . (2) 
a a 

The conjugate momenta of the U's also reside on links and are denoted by E x , i — 1, . . . , N 2 — 1. The following are 
the commutation relations of this set of operators, 

[KAu x ) ab ] = {xu x ) ab , (3) 

[Ei,El] = if* El (4) 

Here A 1 are matrices that represent the traceless generators of SU(N) in the fundamental representation, and / ufc 
are the structure constants of SU(N). We choose a normalization where 

tr (A l A J ) = (5) 



N A -1 



E X abKd = ^UadScb- JjSabScdj ■ (6) 



1 We note here that, in contrast to what happens in the lattice strong-coupling expansion [24| , the baryons are not tied to lattice sites in 
our approach and their wave-functions are determined dynamically. 

2 As discussed below, this will give rise to extra zero modes (in addition to the ones already mentioned above) but because the real ground 
state has vanishing Polyakov loops and no extra zero modes, we do not investigate this issue further. 



3 



Let us now discuss the lattice Hilbert space. A general state |0) is the direct product on all lattice sites, 

|n) = \Q) G ® |n>*. (7) 

Here the first factor is the projection of the state |0) to the gauge field sector, while the second factor describes the 
fermionic sector. Any state is the following direct product 

i^*=n(i fi )*)x- (8) 

Concentrate on the Hilbert space of each site : the "lowest" state is the no-quantum "drained" state |dr) a , defined as 

C a |dr) x = 0. (9) 

Applications of various ipx" create the corresponding quarks on that site, 

\aa) x * = ipt aa \dr) x . (10) 

For the free theory, a = 1 correspond to creation of positive energy excitations and a = 2 corresponds to creation of 
negative energy excitations. To use the usual quark-anti-quark language we write for each site and color 

b' creates a quark and eft an anti-quark. From Eq. we see that b and eft annihilate the drained state, which means 
that this state is empty of quarks, and filled with anti-quarks. The local baryon density operator is 

iv 



N ' lx N 

a— 1 a 

According to Eq. (| 1 2[> . the baryon number of the drained state is —1, corresponding to filling the site with anti-quarks. 
The vacuum |0) is the state with no quarks and no anti-quarks. This state is the filled Dirac sea on a single site and 
obeys, 

6|0) = d\0) = 0. (13) 

The baryon number of this state is B = 0. Because of the Pauli exclusion principle we cannot put too many fermions 
on a single site. For a single flavor theory the maximum number of local baryon number will be 1, and is found only 
in the state | filled) 

6 t |filled) = d|fiUed) = 0. (14) 

Below we will see that gauge invariance puts more restrictions on the single site Hilbert space in order that it be color 
neutral. 

Moving to the gauge Hilbert space, we also write it as a direct product of the form 



P) G = j\m G ) x . as) 



Next, we denote the state with no electric field E by \0)g> 

E l \0) G = 0, V*. (16) 

Any application of the link operators (U x ) ab on \Q)g, creates states which correspond to flux lines on the link x. The 
state |0)g is the only state with no flux at all. Using the fact that the electric field operators generate a SU(N) 
algebra, one can distinguish between the different quantum states created by the link operators as follows. Define the 
quadratic Casimir operator 

JV 2 -1 
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It is clear that the flux-less state \0)g is an eigenstate of this Casimir, with zero eigenvalue. Next, the commutations 
of the Casimir with the link operators are verified from Eq. ([3]) to be 



El(U x ) ab =C F (U x ) ab . 



(18) 



where Cp — (N 2 — 1)/2N is the Casimir operator in the fundamental representation. This means that the state 
[(U x ) ab \0)G] is also an eigenstate of (fT7|) . with eigenvalue equal to Cf- One can now classify the states in \Q)g 
according to their E 2 eigenvalue. The result is a Hilbert space with a ladder- like structure. The lowest state is |0)g, 
with zero flux, and is a singlet of SU{N). Repeated applications of the gauge field operators (U x ) ab create states with 
higher and higher values of flux and the operators E 2 measures the flux on the link x. Indeed we shall see shortly 
that it is proportional to the (kinetic) energy of the gauge fields. 

To complete the picture we now discuss gauge invariance. First recall that the starting point of this formalism was 
to choose the time-like gauge. This leaves only time-independent gauge transformations as a symmetry. The fermion 
operators belong to the fundamental representation of the gauge symmetry and transform as 



with V € SU (N) given in general by 



V x = exp 



N-l 



E ^ 



(19) 



(20) 



Using the anti-commutation relations |T]), one can show that the quantum operator Vf that realizes Eq. (fT^) in Hilbert 
space as 



(21) 



is 



exp 



Af 2 -1 



E PA^M 



x i=X 



The gauge fields transform according to 

(U x ) ab - V G (U x ) ab v G = (v x u x v^ +1 ) 

\ / ab 



(22) 



(23) 



Using the commutation relations in Eq. f3j, one shows that the quantum operator Vg that generates these rotations 
is given by 3 



Vg = exp 



N 2 -l 



x i—1 



EE^N-(^U^ 1 



(24) 



where here the matrix of operators (E/^j[' ) is the link matrix in the adjoint representation, i. 



([/^■V. =2tr (XU x XUi). 



(25) 



3 Here we take 9 X to be periodic in x. 

4 Some useful properties of the operator 
and that J2 k (^ Adj ) - fe (^ Adj ) , fc = Sy. The latter relation means that for the operator identity (C7 Ad j- 1/tAdj.j 



(U Ad i^ are that its representation in Hilbert space is Hermitian, (t/ Ad -i') = ^(t/ Ad -i ) J , 



need to define ([/tAdj.j = (f7 Ad J-) 
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Putting Eq. (j22f . and Eq. ([24]) together, we see that the operator that induces gauge transformations is 



V = cxp 



(26) 



where p l Fx , and are the color charge densities of the fermions and of the gauge fields. These two quantities are 
given by 

Pfx = ri a K b ipl (27) 

pb x = K-(u^{) El_ x =DEl (28) 

Since the lattice Hamiltonian is gauge invariant, we know that the generators of the gauge transformations commute 
with the Hamiltonian 

[H,p Gx -p Fx ] =0, Vz Vs. (29) 

This means that we can choose to work with a basis that block diagonalizes p l G — p F . This breaks the Hilbert space 
to separate sectors classified by their eigenvalue paternal, x = Pg ~ Pf 011 the lattice. Each set of values /^external x 
describes a different physical case, with a different external distribution of color charge (that can correspond, for 
example, to infinitely heavy quarks etc.). To describe the physics of zero external gauge fields, we work with the 
choice Pextemai = 0- Working in this subspace means that all physical states must be color singlets, since all gauge 
transformations are trivial. This means that in this sector the following equations hold as operators identities 

(' ; v: r) 4-i = pf x - m 

Finally we write the lattice Hamiltonian H of the 1 + 1-dimensional SU (N) gauge theory with one flavor of fermions. 
H is given by (in one dimension there is no magnetic field and so the plaquette term is identically zero) 

H = H E +H F . (31) 

Here He is the electric term, a sum over links x of the SU (N) Casimir operator 

H G = £jr(El) 2 . (32) 

x=l 



Next is the fermion Hamiltonian, 

L 

2 



H F = -\Y,^ aa (^) a pU x , ab ^ b f +1 + h.c. + mJ2^ aa Mapflf- (33) 



Here we choose a particular representation of the one-dimensional Dirac matrices using the Pauli matrices and 

periodic boundary conditions on the fermions, i.e. ipx = V'x+Ls > an d "01 = V^+l ■ 

Superficially, the first term in the Hamiltonian is symmetric under a U{1)r x U{1)l chiral symmetry which is 
explicitly broken by the mass term to the vector U(l). One can, however, spin diagonalize the fermions by writing 
ip° -> {e- llxa ^l 2 al) a ipg and see that 

HF = ^Y,^ aaU ^l a +1 +hM. + m Y / hi) x ^ x aa (<r 3 ) a ^f- (34) 

X X 

Here, while the first term is invariant under a U(2) group, the mass term is invariant only under a U(l) x U(l). In 
our actual calculations we will dro p th e second component of this spin-diagonalized Hamiltonian (the a = 2), and 
thus work with staggered fermions [231 ■ To go back to Dirac fermions is easy and for our purpose it is useful to note 
that the continuum chiral condensate and baryon number local densities, at position X, Oipip) (X) and (jp y) (^0> 
are given, up to overall renormalization factors, by [2^ | 

tyty) PO = 4 (l&k + 1>t-i1>*-i) , (35) 
(W) (X) = ±= (yl^ - ^_^ x _i) , (36) 
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where we take X/a = x, and x denotes an even site. 5 



III. AXIAL GAUGE FIXING 

In this section we show how to fix the axial gauge. Since axial gauge fixing in the Hamiltonian approach is less 
familiar than it is in the Euclidean approach we begin with an explanation of the general strategy. 

Temporal gauge fixing left us with a Hamiltonian that is invariant under time-independent gauge transformations 
discussed in the previous section. With no external charges the generator of such transforms needs to vanish on 
physical states and this gives rise to a set of local and global Gauss-law constraints that the quantum fields and their 
conjugate momenta need to obey. Of these constraints, the local ones can be solved (and will be solved in the next 
section) and consequently a large subset of the gauge fields' conjugate momenta is written in terms of the fermion 
color charges. The fact that these momenta become non-dynamical implies that their conjugate gauge fields are not 
physical and can be removed from the Hamiltonian. These are the fields that, in the action formalism, can be 'gauged 
away'. 

Indeed, from the path integral formalism we know that in a compact system of one spatial dimension almost 
all gauge fields can be gauged away. The only gluonic modes that remain are those corresponding to a constant 
spatial gauge field. Furthermore, one can gauge away all but the N — 1 independent eigenvalues of this zero mode. 
Anticipating a similar scenario in the Hamiltonian approach, one expects that the following fermionic Hamiltonian 
will be the remnant of the Hamiltonian Eq. (|33[) after the axial gauge fixing : 

H> F = -iJ2 V4V 3 ^ C+i + h.c. + mJ2 V£Vi C. (37) 

X X 

Here we denote by e liPa the a th eigenvalues of the constant spatial gauge field operator. Clearly, to get Eq. |37j) from 
Eq. (|33|) one needs to remove all but the eigenvalues of the spatial Polyakov loops operator from the system. Using a 
sequence of a change of variables, this is easy to do in the path integral approach. But in the Hamiltonian there seems 
to be a conceptual difficulty with such a process : the degrees of freedom we wish to gauge away are represented by 
quantum operators and it is not clear how to 'gauge them away'. 

To proceed we choose to generalize, to the lattice, the formalism constructed in the seminal Ref. [23| , The general 
idea is to find a Hilbert space realization of a unitary operator T that will rotate quantum states into a basis where 
the Hamiltonian looks like Eq. (|37| . To do so we follow the path integral picture as a guide : first we define the 
spatial Polyakov loop operator P to be 

N 

Pab = ^2 U\,acU2,cd- ' 'Uh^zb , (38) 
c,d,- ■ ■ .z—1 

and second we define the 'eigenvalue operators' <p a through 

Pab=Y, S U P ) etLsMP)S ^(P)- (39) 

c 

Here the matrix of operators (S(P)) ab is a functional of the operators P ab and is defined implicitly through Eq. (|3D[) . 
Next we ask to find a form of T that will induce the following transformations 

*r x ^ = J2 v ^ b ^l (40) 

b 

FTPl a ^=J2^ bV l b a, (41) 

b 

U' xab =J 7 U x ,ab& =U x ,ab, (42) 



< - 


- TP*" 


- 




Ux.ab ~ 


- K„ 


u x,ab 





't _ fj]\ f \ _ T]\ 
b ~ •> u x.ab- r ~ u x,ab- 



(43) 



5 Here all fields arc still dimcnsionlcss. The continuum expectation values ^ 1,70 ^L— are given by dividing the r.h.s. of Eqs. (J35J 

(136 I t by a-v/ g 2 N (g has dimensions of mass in 1 + 1). 
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Crucially, here the operator (V x ) ab is chosen to be the following functional of the gauge fields operators (U x ) ab , 

N 

V x , ab = £ U Uac Ul_ 2M . . . u{ yz Sl b (P) *. (44) 



cde---yz— 1 



This guarantees that V obeys 



E(^) ob (U x ) bc (V x+ i) cd = e i ^5 ad 



(45) 



be 



as an operator identity. The form of T that induces Eqs. (|40|) (|42|l is then written as 



(46) 



where G is again an operator in Hilbert space that depends on the gauge field operators U x _ a t and that is defined 
through its following relation to the operator V : 



V x _ab = 



exp ire* ({[/}) 



(47) 



- ab 



Note that T is not a gauge transformation (i.e. its not of the form Eq. ([26)) ). 

The end result of Eqs. (|4"0|) -([47 |) is that when we rotate the Hilbert space by T or, equivalently, conjugate the 
fermionic part of the Hamiltonian by T we indeed get Eq. (137|) : 



H F -> H' F = TH F T^ = -i^^v 3 e llp -ip a x+l + h.c. + m^^l a a l ip a x . 



(48) 



This change of quantum basis is how the process of axial gauge fixing works in the Hamiltonian formalism. 

We now proceed to find the transformed version of He, of the commutation relations, and of the Gauss-law 
constraint. We begin by defining the transformed electric fields E' via 



Ej ee TE X T\ 

which allows us to write the transformed electric field Hamiltonian as 

,2 



H E - H' E =TH E P = y £ (4*)' 



(49) 



(50) 



It is clear that the commutation relation between the E' fields and the U fields are the same as in Eq. ([3]) and Eq. ((4]) 
(this is so because of Eq. (|42p). Next, we turn to transform the Gauss law constraint. For that we conjugate Eq. (150)) 
by T 



This equation, in terms of E*, reads 

e'J - (u x -i) e'i, = f Pf>x ^ = J2 y ab t = £ Kb € a ' € vL a v x 



bb' 



(51) 



(52) 



Multiplying this equation by \ l cd , summing over i, and writing it in terms of the adjoint representation of the operator 
V, we get 



(53) 
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Finally, if we define 

K" = (V x Adl ) lk E'J, (54) 

and 

E'J = (V x Adi -) jk E x \ (55) 

then we get 

^-(^dK-i).. /'••', />',.,• (56) 

Using Eq. (|44|) this gives 

<M e ~^ dj ^i = PV (57) 

Eq. ([57| is the starting point to the discussions in Sections HVtfvl : first, in Section Hvl we solve Eq. (|57|) and write 
the fields E x l in terms of the color charge densities p Fx - This cannot be done to all of the components of E , and a 
subset of these remains independent of the fermionic charge densities. The way we treat these remaining components 
is explained in Section fVl 



IV. RESOLUTION OF THE GAUSS LAW CONSTRAINT 



In this section we resolve the Gauss law constraint of Eq. (|57| and write the electric field operators E in terms 
of the fermion color charge densities pf- For that purpose we first specify the basis of the color group generators A\ 
We choose the N — 1 traceless generators < \ 1=1 > 2 >--->JV-1 to span the traceless Cartan sub-algebra, and the remaining 
N(N — 1) generators to have no diagonal entries. In this basis we see that 



Sij if i,j e [1,7V - 1], g 

if i € [1, iV— 1] but notj, and vice versa ' ^ ' 



As a result, if we focus on the Cartan sub-algebra (let us denote its generators by i = I), and Fourier transform the 
operators E" and p according to 

El 1 = £ e^E'; 1 , (59) 



iU (60) 



v p 



with p — '^Yp- ;n = 0,l,...,L s — 1, then we can resolve the nonzero momentum components of E L : 

K 1 = for (61) 

For p — Gauss law cannot be resolved, and instead becomes a constraint on the fermionic global color charges : 

V/=l,2,...,iV-l : Q 1 = J2^i Xl ^ = °- (62) 



The N — 1 requirements of Eq. (|62|) can be written as 

Va = 1, . . . , N : £ V4°V£ = independent of a, (63) 

X 

and since the baryon number B is equal to -4 J2 X a V4 a V>x ~ L s (see Eq. l[T2"j) ) then Eq. becomes 

^Vi a C = S + i s , Va= 1,2,...,JV. (64) 
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For the rest of the generators i ^ / we can indeed resolve Gauss law even for zero momentum : multiplying Eq. (|57|) 
b y Kb witn a ^ b, and defining E xab = Y^tyi Kb E x\ and PF.x.ab = J2i^i KbPF,x' one nnds tna t the Fourier 
components of E x a ^ b obey 

PF,p,ab 



E 



p,ab 



1 — e -i(p+<Pa-¥b) ' 



(65) 



(here the Fourier transformations of E x a _ ib and p x . a ^b are defined in a way similar to Eqs. (|59l - |60"f ) Note that we 
work in a Schroedinger picture where ip a is a c-number. Also, in Eq. (|65|) we assume the absence of states with 
fa — fb + P = 0. Since p is quantized in units of 2it/L s , this assumption will indeed become true - see the discussion 
below on the importance of the Jacobian of the curvilinear coordinates (p. 

We can now use the relations between E and pp , and write the electric Hamiltonian in terms of the fermions and 
the gluonic zero modes. For that we use the hermiticity of E' x and get 



H' = 



9' 



E«f yX'^v 9 ^E E * iiE * ^Xv t v 



/ 



f£_ 

2 



2 x ^ Px.ab Py,ba 



i{y-x)p 



L s ^ 4sin 2 (H±M; 



N-l 

i=i 



(66) 



Here, by the primed sum, we mean that one must sum over all a, b and p for which the denominator is nonzero. This 
restriction is a direct consequence of one's inability to resolve all of the Gauss law constraints and will become the 
source of the principle value prescription often used in the 't Hooft model. 

V. REALIZATION OF THE ZERO MODES IN THE SCHROEDINGER PICTURE 



The resolution of Gauss law left us in the glue sector with a restricted set of dynamical degrees of freedom the 
N — 1 eigenvalues of the spatial Polyakov loops, ip a , and their N — 1 conjugate momenta E^ =(j . In this section we 
focus on these operators to which we collectively refer as the zero modes. Specifically, we calculate their commutation 
relations and also find a simple way to realize Ep =0 in the Schroedinger picture, where ip a are c-numbers. For that 

purpose we first calculate the commutation relations of E x % and the operator 



(Py) ab = (UiU 2 ---U y ) ab 



(67) 



The result is easy to obtain and can be written in terms of the operators S and (p (see Eq. (|39|)). and the operator V 
(Eq. (gl), 



E'l (P„)J = 6{y - x) e **l2v V^'SPy)^ 



(68) 



Here 6{x) = 1 if x > and zero otherwise. Next we use the definition of E 1 in Eq. and after some algebra find 
that (here note that we focus only on the Cartan sub-algebra) 

E7, (Py) a b} = 0{V - s) K SP v ) ab . (69) 

Since the dependence on the site index x in the r.h.s. of Eq. (|69p is trivial, we can easily write the commutation 
relation of the zero mode E p L — -j= J2 x E,J , and we find 



Ppio' (Py)ab 



In particular, for y = L s , we use Eq. (155)) and get 



EpLo-, Pab 



(S* X 1 SP y ) ab 



If we now define £ a = J2i Ka E P lo> we § et 



[E c ,(tfe^S) ab ] = l^/L s \st c e^>S cb -±P a 



(70) 



(71) 



(72) 
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In Appendix [X] we show that 



£a = £l (73) 

[£ a ,£b] = 0. (74) 

and we note in passing also that £ a = 0. 

The form of Eqs. ([72|) -(|74 |) leads us to write the following realization of the £ operators in the Schrocdinger picture 



S logA 2 (<^) 1 ^ Slog A 



5tp a N Sip, 



(75) 



N Sf c J 2 
where here A 2 (tp) is the well known Vandermond determinant 

A 2 (M) = J] ^ 2 (l.^=^ ) . (7<i, 

a<b 



It is easy to check that the first term in Eq. ([75]) indeed satisfies Eq. (|72p : the (p a =i,...,N behave like the coordinates of 

N particles, while ^ jj J2b behave like their momenta, relative to the motion of the center of mass coordinate 

ipcm. = ^2 a if a - This separation into relative coordinates and a center of mass is anticipated in the SU(N) case where 
ipc.m. is not a true degree of freedom. 

It is also not hard to understand the origin of the second term in Eq. (|75|) : while it trivially obeys Eq. (|72|) . it is 
necessary for Eq. (|73p . To see this recall that in the Schroedinger picture the measure of the ip a coordinates is not 
flat - it is given by the Haar measure dP over the spatial Polyakov loop P (which is the only remnant of the gauge 
field degrees of freedom). Since the Hamiltonian only depends on the eigenvalues ip a L s of P, then we can replace 



J (jl^ A 2 {{tp})5{tp c . m ). 



dP^ I \\\cbp a ] A z ({<p})S(<p c . m .). (77) 

SU(N) 

Because of the A 2 factor in the measure, the simple derivative operator — ij- will not be hermitian (its action to the 

right will differ from its action to the left by a derivative of A 2 ). It is easy to check the particular choice of in Eq. ((75|) 
fixes this and makes Eq. (|75|) hermitian. 6 

Our final step will be to simplify Eq. (|75|) even further. For that we define a new wave function f llcw for the 
curvilinear coordinates ip that will make their measure flat. Specifically, we write 

In terms of the new wave function the Vandermond disappears from the measure and the kinetic term becomes the 
simple quadratic form (here we rescaled ip — ► ip/L s ): 

This will be our final simplified form for the zero mode contribution to H E . 



VI. RECAP : THE HAMILTONIAN AND THE RESTRICTIONS ON THE HILBERT SPACE. 

To conclude let us write the lattice Hamiltonian H in terms of the operators that cannot be gauged away. It 
describes the interactions between the lattice fermions ip and the eigenvalues of the spatial Polyakov loop <p, and is 



The integral over <p can be either restricted to obey ip a > fb for a > b, or can be simply unrestricted. The difference between these two 
choices will be reflected by whether the operators S(P) include permutations of the eigenvalues or not. 
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given by 

H 



Pf,: 

Here the primed sum in Eq. ([55]) means we that we should not sum over terms whose denominator is zero (for the 
derivation of this restriction see section HV)l . From here on we discard the lower component of the fermions, and so 
discuss the staggered fermions formulation of lattice QCD [27| (this also means that we need to replace the r.h.s. of 
Eq. (|64|) by B + L s /2, and choose the number of lattice sites, L s , to be even). The reason we choose to work with 
staggered fermions is two-fold. Firstly, these one-component Grassmann variables are simpler to treat in our formalism 
and the computational effort involved in the numerical minimization of their classical Hamiltonian (for details see 
[Til) is a factor of two smaller than for the naive prescription of fermions considered above. Secondly, unlike the 
situation in four dimensional euclidean calculations, staggered fermions in the one dimensional Hamiltonian approach 
are 'un-doubled' and correspond to a single Dirac fermion in the continuum limit. 

Finally, the Hilbert space of the system obeys the following : 
• The operators J2 X ^x^x are > f° r eacn vaiu e of the color index a = 1, . . . , N, constants and equal to B + L s /2. 



Hg + Hk + He, 
g 2 L s ^ { 5 I ^ ,) 



E 

d=l 



1 N 

-Y 



Sfd N Sip, 



€ a e Wis r x+ i + h.c. + m j2 (-if 4 a r x , 



Pf.x Pfa 



1 ( ( Va - Vb )/L s +p\ ' 



4sin^ 



(80) 
(81) 

(82) 
(83) 

(84) 



VII. DIAGONALIZATION OF THE HAMILTONIAN : THE COHERENT STATES APPROACH 

In this section we use the coherent state approach to study the ground state properties of the Hamiltonian in its 
large- ./V limit. This approach is summarized in Refs. [25| and we refer to these papers for detailed discussions, while 
here we only describe its strategy. 

The main paradigm that underlies the coherent state approach is that QCD, in its large-iV limit, becomes a 
classical theory. 7 This means that instead of diagonalizing its quantum Hamiltonian, one can instead minimize a 
corresponding classical function TL, referred to as the 'classical Hamiltonian'. The minimization is done with respect 
to a set of coordinates, C, that corresponds to expectation values of gauge invariant operators. The coherent state 
approach provides the mathematical prescription for calculating the function TL{C) and it is shown to be given by 

H{C) = (C\H\C). (85) 

Here \C) is the so called 'coherent state' obtained by applying a unitary color singlet operator, that we denote by 
U{C), on an arbitrary color singlet element |0) of Hilbert space, the so called 'reference' state : 

\C)=U(C)\0). (86) 

The operator 11 is a functional of gauge invariant operators like spatial Wilson loops of contour T 

tr W r , (87) 



7 Here we implicitly assume that we take the large-JV limit when all other parameters are fixed. These includes the volume L, the 
temperature T, the mass m q , the 't Hooft coupling g 2 N, and the baryon number B. 
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as well as spatial Wilson loops decorated by a single insertion of the electric field matrix operator (E x ) ab = E % A a fc, 

i=i 

at point x along the loop, i.e. 

tri^grWr, and trW r £*er. (88) 



(the two operators that appear in Eq. are different since the electric fields do not commute with the gauge fields). 
Also, in the presence of fermions, IA depends on string-like operators of the form 



(89) 



where here U x ^ y is a string of gauge field operators connecting site x with site y. The generic way U(C) depends on 
these operators is 



U{C) ee U F xU G 
Uf = exp 



exp 



xy 



r xer 



(90) 
(91) 

(92) 



form a closed algebra, the structure of Eqs. 



where C F V ,C G , and C G r 1 2 are c-numbers. 

Since all the three type of operators defined in Eqs. 
guarantees that U(C) is an element of a Lie group referred to as the coherence group. Indeed C F V ,C g, and C G \ 2 
parameterize this Lie group and furnish coordinates on its manifold. More precisely, the fermionic part of the coherence 
group is parameterized by the fermionic coordinates Cf , and the gluonic part by C G and C G ,i.2- s 

The values C m in that minimize Ti.(C) then determine the values of all gauge invariant observables in the ground 
state of H, so for example, the ground state energy £ g . s , is given by 



N 



Jim Eg. B . = H(C = C mi „), 



(93) 



etc. 



In the large- N limit of QCD with the number of flavors, Nf, kept fixed, the minimization of TL process proceeds in 
two steps : one begins by minimizing the leading 0(N 2 ) contribution of H(C). At leading order the latter is equal to 
the classical function TLq given by the classical Hamiltonian of the pure gauge theory 



H G = {0 g \W G H G U G \0 



(94) 



where the operator H G is given in Eq. (|5T|) and |0 ff ) is a reference state in the pure gauge system. This step sets the 
value of the expectation values of the gluonic color singlet operators such as the spatial Polyakov loops. For brevity, 
let us refer to these expectation values by the generic symbol V, and by Pmin to their value at the minimum of H G . 

Next, one minimizes the O(N) contribution to H. To leading order it is given by calculating the expectation value 
of Hk + H G from Eqs. (|52"|) and (|53"|) in the subset of coherent states whose gluonic coordinates have already been 
determined by the gauge dynamics. This contribution is thus given by 



H F (C F ) ee (0 F \Ul (H K + H c ) U F \V F )\v=v a 



(95) 



with | Of) denoting a fermionic reference state. This two step process reflects the dominance of the gauge field dynamics 
over the fermion dynamics. Indeed the back-reaction of the fermions on the gauge fields comes from terms which we 
do not consider in this work and that are subleading in 1/iV. This is true as long as there is no enhancement of these 
1/N terms by massless modes, which we assume to be the case. 9 



To ensure that U(C) is unitary, the coordinates C need to obey certain conditions and in the next section we make these explicit. 
We thank V. P. Nair for pointing this to us and to L. G. Yaffe for a discussion related to this. 
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For our purposes it will be easier to simply calculate TLf{Cf) and substitute the set of expectation values V by its 
value in the exact ground state of Hq- The latter is known analytically and we repeat its derivation in Section [VII Al 
The solution of the fermion sector, however, is done as described above and we calculate the classical function Hf(Cf) 
in Section ETUI 

We now turn to make the following remark on the way 7i is calculated and minimized. As mentioned above, the 
coordinates C corresponds to expectation values of different color singlet operators. More accurately, one can focus 
ones attentions to expectation values of only 'non-factorizable' operators such as the string operators of Eq. (189p . 
Naively, however, one might expect that TL should also be minimized with respect to coordinates that correspond to 
expectation values of 'factorizable' operators such as 

^ X , ^ fat U^ylPy) X 1>,) , • • ■ ■ (96) 

This, however, is incorrect : at large- N the expectation values of such operators factorize, and for example, 

(V4V* >< ^yi>v) N =°° (V4V>*) X (^ y ). (97) 

As a result, the expectation values of factorizable operators become determined by the expectation values of the 
non-factorizable operators and should not be thought of as independent coordinates that parameterize the coherent 
state manifold or the point within that manifold that represents the ground state. 

The unique role of non-factorizable operators operators is also reflected by the structure of the unitary operator of 
Eqs. (|9"0H9"2"|) which depends only on such operators. Indeed in Ref. [25[ it is shown that to generate the whole coherence 
group, it is sufficient to include in its algebra only the non-factorizable operators that appear in Eqs. (f9"0H9"2")) . Adding 
other operators greatly complicates the algebra of the coherence group and is unnecessary. An important result of the 
discussion above, which we wish to emphasize, is that the coherent state \C) defined in Eq. (|86p is sufficiently general 
to look for all possible large- N ground states. Put differently, the form in Eq. ([86| does not correspond to assuming 
an ansatz for the ground state of the gauge theory. 

In the next two subsections we show how to implement the program outlined above in practice. We begin in 
Section IVlI Al with the treatment of the pure gauge case that allows us to calculate the expectation values of the color 
singlet gluonic operators, namely the traces of different powers of spatial Polyakov loops. The pure gauge Hamiltonian 
in our two-dimensional case is sufficiently simple that we can do so exactly and for any value of N. The next step is 
to apply the coherent state approach to the fermionic part of the Hamiltonian and we do so in Section IVII Bl 



A. Solution of the gauge sector 

The solution of the pure gauge sector is well known [28| . The starting point is to notice that Eq. lj78|) means that 
the gauge wave functions need to be anti-symmetric to an exchange of two angles, and for odd values of N, periodic 
in 2n: 

i &(<Pl,V2,---,<Pk,>--,<PU---,<PN) = -^(fl,f2,---,fl,---,(fk,---,fN)- (98) 

4 f (<pi,<p2,...,<Pk + 27r, ...,<£iv) = +$((pi,ip2,---,ipk,---,<PN)- (99) 

These properties, together with the form of the Hamiltonian, tell us to think of the ip a as the positions of N non- 
relativistic fermions with mass 2/g 2 L s , moving on a circle with periodic boundary conditions and a fixed center of 
mass. The single particle wave functions of such a system are the plane waves e lvn with n = 0, ±1, ±2, . . ., and the 
ground state wave function of this A^-fermion system is, up to a phase, the N x N slater determinant obtained by 
occupying momentum states distributed symmetrically around zero and limited by 'Fermi momenta' uf = (N — l)/2: 

*({*>})= det e ib ^. (100) 

-TiF<(t,K-fUjr 

This wave function Eq. (|100|) is an eigenstate of Hq with eigen-energy that has the large- N limit of L °J£ x g 2 N . 

It is easy to show that this determinant differs from the Vandermond determinant Eq. (|76p by a phase that depends 
only on the 'center of mass' of the fermions, ^ c ip c , and in SU(N) we set this phase to zero. The result is 

|*(M)| 2 = a 2 (M), (ioi) 
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and so when we calculate expectation values of gauge invariant gluonic operators 0({ip}) we need to perform the 
following integral (here we set the normalization such that (1) = 1 and denote ^ c tp c by tp c .m.) 

{6) G = /^A 2 (M)^ c . m .) o(M), (102) 



B. Solution of the fermion sector : large- N coherent states 



We diagonalize the fermion sector in a variational manner. As discussed above, the most general form of the 
coherent state \C F ) is given by 

\C F ) = U(C F )\0) =exp l-i E CF^e^y-'W-W I |0), (104) 
\ xez Ls y ez J 

and parameterized by the infinite dimensional hcrmitian matrix C F . The state |0) is a fermionic reference state which 
we choose to be a state annihilated by V4 for a subset of lattice sites that are full of fermions and that we denote by 
S. More precisely |0) is defined to obey 

4,J0>=0 ; VxeS. (105) 

On the complementary set of sites, x £ S, the reference state is annihilated by %j) x . This means that we choose 
B + L s /2 of the lattice sites to be full of baryons and the rest empty of baryons. This choice is convenient for our 
calculation, but the results are insensitive to it : the only thing that matters is that the overall baryon number that 
|0) contains is B + L s /2 (such that the renormalized baryon number is B). 

We emphasize here that in the exponent of Eq. (|104|) . the index x runs over 1, 2, . . . , L s , while the sum over y is 
unrestricted : y = 0, ±1, ±2, . . . , ±oo. The identification 

r x =r x+kLs ; kez, (106) 

means that Cf needs to obey 

c-- iL:, - iL ■ kez, (107) 

in order to have U(C F )W{C F ) = 1. 

To find what is Cf we minimize the classical Hamiltonian TL defined by 

H F = [ d<p (C F \H K + H C \C F ) A 2 {ip) x 8 (</w.) , (108) 



in the space of all possible choices of C F . The calculation of (C f \Hk + Hk\C f ) is somewhat lengthy and we postpone 
it to Appendix [Bl 

The resulting expression can be written in terms of the following 'density matrices' p%. y and p% y with x, y G [1, L s ] 
and q G Z 

Z+P L S , X ^ tCF y,z+(,+P)L s : (10Q) 
ly, (HO) 

In terms of p% „. and p% ai the classical Hamiltonian is 

(c f \h k + h c \c f ) = Yl [(44+i+ c - c -)+ m (-ir&l><E e< "" 

qez 

+ 1- y y y p *y p y* e + o(i/n). 

ALs Ml. 4sin 2 ((^-^)/L s +p ) /2 



p q 

rxy 


- E ^ Cf ) 










p q 

rxy 


= Sq.oSxy — P 



(111) 



15 



Note that in contrast to the form Eq. (|83|) , here the second term does not contain any contributions from the a = b 
terms. They are shown to be subleading in Appendix [Bl 

To obtain a compact form for the classical Hamiltonian Tic we define Pk to be the Polyakov loop operator that 
winds k times around the torus as 

p ^^E elVafe ' ( 112 ) 



and use the identity 

1 



\Q\e~^ Qk , (H3) 



4sin^((£ : + *e)/2) 2 Q/ 

with e > regularizing the pole of the left hand side, to re-write the second term of Eq. (|111|) as a sum over terms 
that only contain powers of e l ^ Pa ~ Vb '' LB . Indeed, substituting Eq. (|113[) into the second term of Eq. (jTTTJ) , we find 
that the terms contributing to the sum over the momentum variable p = 2nl/ L s can be isolated and read 

^ e - iMix -y- Q)/L , =Ls ^ 5y _ x+Q QLs . (114) 

This allows us to substitute Q in Eq. (|113p by x — y + QL S . Since the dummy summation variable Q obtains all 
possible integer values, we drop the 'bar' from its notation in the rest of the paper. Thus, the final result of the 
manipulations in the last paragraph is that Eq. (|llip gets the following form. 



Hf(p)/N = 



q£Z 

9 2 N 



+ c.c. +m(-l) x pl a 



E E pi v P^\x-y + QL s \e- elx - y+QL ' ] ({P Q - q P- (Q - ql) ) G - ^{P q >- q ) G ) , 
xy£Z Ls Q qq 'ez 



(115) 

where here we define the gluonic exception values (, }g in Eq. (|102|) . 



1. Remarks on the form of the classical Hamiltonian 

Eq. (|115|) above is quite an important ingredient in our paper and so let us now pause here and make the following 
remark on its form and its implications. What Eq. (|115[) tells us is that the fermionic properties of the system, 
which are represented by the density matrices p% y , self interact as well as couple to the Polyakov loops P q . The 
expectation values of the latter are determined by the gauge dynamics and feel no back reaction from the fermions. 
Thus, performing the minimization in Section [VII Al correctly, and determining the properties of the gluonic vacuum 
in a consistent way, is crucial to get the correct fermion dynamics. For example, in the Sections IVIIII and IIXI we 
emphasize that mistreating the glue sector, which is what one does when one ignores the zero modes, leads, through 
the way Eq. ()115|) couples the glue sector to the fermion sector, to erroneous results for various fermionic expectation 
values. 

In particular, one can see that if the expectation values of the different Polyakov loops are incorrectly chosen to be 
unity, 

P q = 1, Vq, (116) 

then by using large- N factorization of the double trace term (PQ- q P_(q_ 5 ')) n = oc (Pg_ g ) (P_( Q_ q ^) , one finds that 
the second term in Eq. (|115p is strongly dependent on L s (this is shown explicitly in Section [IX A[) . This volume 
dependence is contradicting general arguments on large- N gauge theories such as those of Ref. [3j . 



2. Further manipulations of the classical Hamiltonian and preparing for its minimization 



In this section we make further important simplifications of Eq. (|115[) that will also allow us to minimize it (see for 
example Ref. [l6| and the next sections). 
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For a U(N) theory one drops the delta functions in Eq. (|102[) . and the gluonic expectation values (Pq)c and 
(P q P q ')c are calculated in Appendix [B] (see Eqs. IB16IIB17)) . In the SU(N) case the integral in Eq. (|102[) changes 
only by restricting the sum ^ ip a to be zero, and so we expect the Polyakov loops to be the same in the large- N limit 
(we show this explicitly for (P 9 )g in Appendix IB"]) . Using Eqs. (|B16|) - (|B17|) we see that, as expected, the leading 
contribution to Eq. |TT5|) is from the q = term in the first line and from the q = q' = Q in the first term of the 
second line. This gives us 10 



H F (p)/N = ]T (™pk+i + c.c)+m(-l) B /4 



xez L 
9 2 N 



E E - v + QL.\< " V+QL ■ 

xy£Z Ls QeZ 



(117) 



Since the index Q runs over all integers we can write 

r7r dp 



P% = I ^P X yiP)e ipQ , (H8) 



which, re-using the identity Eq. (|113[) gives the form 



H,.(//)/.V / J E {(-^.x+iOO+cc.) + ml -!)>,,(/-) 



«? 2 A f f dp dp' 1 yv ^ p xy (p)Pyx(p')e i2 ^ x - y V L 
xy€Z Ls I 



Up Up 1 v - x - 

2^2^ZT ^ ^4sin 2 (0-p')/i s + 27rZ/i s )/2 



(119) 

The pole at p — p' and I = L s in Eq. (|119[) might seem alarming and tracking back its source to Eq. (jllip one 
finds that it is the double sum over the color indices a and b that appears in the second term there. Specifically, 
terms in that sum for which \a — b\ is small are causing this divergence. For these terms the argument of the sine 
in the denominator can be small; for example if p = 2n and the Polyakov loops have zero expectation value in the 
gluonic ground state, then heuristically ip a ~ 2-na/N. (This is not a gauge invariant statement, but for the current 
discussion this subtlety is not important. In the rest of this paper we take great care to avoid such statements when 
it is important to do so.) Thus, we see that when \a — b\ < 0(1), then the argument of the sine in Eq. (jllip is very 
small and this corresponds to the pole in Eq. (|119[) above. 

This pole is the source of the IR divergence in the usual treatments of the 't Hooft model, and that is usually 
resolved with the ad hoc principle value prescription. In our case this is not neede. The a = b terms were excluded 
from the sum in Eq. (|111[) . The way this restriction emerges in Eq. (|119p is through certain conditions obeyed by the 
density matrices p® y and p® y which we present in Appendix [B] (see Eqs. (|B10[) — (|Blljl ) and that reflect the unitarity 
of the operator U{Cf)- In the language of p X y{p) these conditions read 

E P X y(p)Pyz(p) = 0. (120) 

yez Ls 

Thus the divergence in Eq. (|119p is removed and we get the principle value prescription. To show that the near vicinity 
of p = p' (corresponding to \a — b\ small but nonzero) is not causing any lower divergences we need to assume a form 
for the p dependence of p xy {p). Instead we have confirmed this numerically whilst minimizing Ti(p) with respect to 

pM- 

We conclude this section by writing Eq. (] 1 1 9[) in a way which is convenient for its minimization. We first solve the 
constraints Eqs. (|B10[> (|Blip on p® y by writing 

M B+LJ2 

p%=mY1 E CW^CW, (121) 

a— 1 n— 1 



10 According to Eqs. l)B16| l- l|B17[l the terms with q = q' ^ Q are also nonzero, but using the representation in Eqs. I|121ll — (11221) we verified 
that they are subleading in N if we take M — > co with or after we take N — > oo. 
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if; f; cwe^cw, (122) 

a=l n=B+L B /2+l 



where for each a = 1, . . . , M, the single particle wave functions </>™(x) ; ra € [1, £ s ] span an L s x L s dimensional space, 

£ C(x)CW=in. (123) 

Eq. (|12ip is nothing but a discretized way to write the most general expression for p and p, that also obeys Eqs. (|B10[) — 
(|Blip and we show that this is correct in Appendix [C] Note that the full space of solutions for p and p is accessible 
only if M = oo. 

In terms of (f>a(x) we get 

W F (0)/7V = i^ £ + +m(-l)X,} 

5 2 iV 1 ' _ ^^e-'t—XftCS'+O) 



4 L S M 



M XV t£ La 4 sin 2 (^(^ + 0/2 



g>N(B + L a /2) 1 ^ 1 

4 ^ 2 ^4sin 2 (1^ + 0/2)' 



B+L.,/2 

= 2 (125) 

n=l 

Here, by the prime on the sums we mean that the terms with a = b and / = L s are excluded. 11 
We now perform the variation of Ti with respect to the M functions <f> a 

S 



7 « " E e ™ E = 0, (126) 



m,b x 



and we find that they must obey the following M coupled nonlinear differential equations (here we use the Lagrange 
multiplier e° to enforce Eq. (|123[) ) 

E K v €(y) = «(*)> (127) 

yez Ls 

with 

K y = + l -(6 y , x+1 ~6 y ^ 1 )+m(-iy6 xy -g 2 Nv a xy , (128) 



1 — , (B+L s /2 \ 

b \ m=i y 



(129) 



and 



1 t-^ e ^ iTTT+'j 

K ab (y,x) = — V ^ -r-. (130) 

;ez is 4sm ^ ^ ml. + T7 y 



11 In our numerical studies we find that, while irrelevant for the minimization of Hf, the last term in Eq. I|124| l is crucial to include in 
order to get the right baryon mass. 
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Since K a b explicitly depends on cf> then the solution of Eq. (|126[) is a self-consistent process. 

Within the space of all functions (j> that obey Eq. (|127|) . the correct solution is the one that has the lowest value of 
Tt. The latter is not equal to J2an € ni since this will count the Coulomb interaction twice. Instead we find 

B+L s /2 M / \ 
Kolution/iV = — E K + E P m + 1)) +m (-1)^ t%{x)4F(x)] 



g 2 N{B + L s /2) 1 1 

^4sin 2 (^(^ + 0/2 



- ■ 2r^,n.- h , ,w„v (131) 



VIII. COMPARISON TO OTHER RELEVANT WORKS 



In this section we wish to discuss the way that our resulting equations differ from those appearing in other works 
that also regularize the 't Hooft model on a finite circle L. 

1. Ref. 22] : this work looked at QCDi + i with the light-cone Hamiltonian, but neglected the curvilinear character 
of the variables (p. Thus the ground state with (p a = was chosen, that in fact has zero measure within the 
correct solution. This is equivalent to assuming the Polyakov loops are all nonzero, and in Section HXl we discuss 
the consequences of such a choice, but in essence it violates volume reduction. In fact, Ref. [22| showed that this 
erroneous ground state leads to a phase transition as a function of L, which clearly contradicts large- N volume 
independence. 

2. Ref. [23| looked at the Hamiltonian of QCD in axial gauge (in any number of dimensions), and in the continuum. 
It poin ted out to the error made in [12] but did not discuss the consequences of this on the 1 + 1 solution presented 
in [22j . In our paper here we generalize the theoretical framework developed in [23| to the lattice regularization. 

3. Ref. 15] : This interesting paper formed one of the motivations for our study. Here the authors showed how 
the phase transition found in [22j disappears when one chooses an appropriate ansatz for the ground state. The 
way this choice is made in Ref. [l5| . however, is not manifestly gauge invariant and breaks a residual gauge 
symmetry of the Hamiltonian. This has the disadvantage of making it hard to argue that the ansatz used is 
exact at large- N, and to construct a gauge invariant fermionic ground state. We avoid this issue in our paper 
by working with manifestly gauge invariant operators and states. This is most naturally done with the coherent 
state approach. Other differences between our paper and [15| is the fact that we use the lattice regularization, 
we do not restrict to zero baryon number or to a particular ansatz for the ground state, and we analyze the role 
of translation symmetry for volume reduction. In particular, in the following section, we show how a 'soft' form 
of large- N volume independence works if translation invariance breaks to one of its subgroups. Finally we make 
the connection with the Eguchi-Kawai volume independence manifest. 

4. Ref. [l|| : To our knowledge this paper is the only one that solves the 't Hooft equation at nonzero baryon 
number for a general ratio \J g 2 N /m. 12 Unfortunately, the authors restrict to study a single baryon only, and, 
like Ref. [22], ignore the curvilinear nature of the <p coordinates. Consequently, at short lengths (or large baryon 
densities) their approach would fail, and exhibit the same phase transition seen in [22| . Clearly this calls for 
revisiting of the topic which we aim to perform in Ref. [Tr| . 

The list we give above demonstrates the usefulness of our current paper : Firstly, it uses the coherent state 
approach, which is manifestly gauge invariant throughout. Secondly, it generalize all former studies to the lattice 
regularization, and extends the study of volume independence to systems with nonzero baryon density and partial 
translation invariance. Thirdly, it opens the way to study the 't Hooft model for arbitrary values of the quark mass, 
spatial volume, and Baryon number. To our knowledge this is the first time these steps are taken. 



Another paper that discusses nonzero baryon number in the 't Hooft model is |xsj | . but there the authors restrict to translation invariant 
dense systems, which as they show, is inconsistent. For translation non-invariant states, the authors discuss only the vicinity of the chiral 
limit. As we already mentioned, the fact that the baryons are massless in this limit makes them very different from the four-dimensional 
QCD case, and it is important to study the m^O case. 



19 



IX. LARGE- iV VOLUME INDEPENDENCE 



It is generally expected of large-TV QCD in 1 + 1 dimensions to be independent of its volume. This equivalence of 
large- TV gauge theories with different volumes was first suggested in [|[ and caused great excitement since it was seen 
to be a potentially easy way to solve large- N QCD on the lattice. Shortly after the papers [Hill showed that this 
equivalence breaks down in the continuum limit for three or more space-time dimensions. This breaking of reduction 
is signaled by the fact that Polyakov loops that wrap around the volume acquire nonzero expectation values. 

It is useful to put this large- TV equivalence in the more general context of orbifold projections between mother and 
daughter gauge theories : in our case the mother theory is large volume QCD while the daughter theory is small 
volume QCD. These projections are expected to become equivalences when the rank of the gauge group, N, becomes 
large. For a review on this topic we refer to Ref. As shown there, these equivalences hold only between certain 
sectors of the mother and daughter theories which are defined to be neutral under certain symmetries. In the original 
Eguchi-Kawai paper it was stressed explicitly that the ground state in both theories needs to be symmetric under the 
center of the gauge group. For SU(N) gauge theories this means that the global Zjy subgroups of the local SU{N), 
that correspond to multiplying Polyakov loops in different directions by a Zn phase e 27rl ^ N , must be unbroken. Thus 
these Polyakov loops are the order parameters of these symmetries and must have vanishing expectation values for 
the equivalence to hold. 

Another symmetry that the ground states of the mother and daughter theories should respect in the volume 
projection case is translation symmetry. This is clear intuitively - how can we describe a theory which breaks 
translations and that as a result has operators with expectation values that depends on the space-time coordinate, by 
a theory that has no volume ? The requirement of intact translation symmetry can of course be anticipated from the 
construction of Ref. [8| , since this is one of the symmetries that define that neutral sectors of a large-to-small volume 
mapping. Indeed this was already explicitly pointed out in the first paper in (25l |. When this symmetry breaks it is 
no longer true that the physical observables in the large volume theory have a one-to-one mapping to observables in 
the zero volume theory. If the attempt to map the large and small volume theories fails, then clearly they cannot be 
large- N equivalent. Unfortunately, we believe that the role of translation symmetry in the Eguchi-Kawai equivalence 
is not fully stressed in some of the relevant literature. The reason, however, is obvious : the QCD vacuum respects 
translation symmetry ! 

But there is at least one physical scenario where one can expect to get broken translation symmetry in QCD and 
that is at nonzero baryon number or chemical potential, where crystals of different sort can form (for relevant literature 
on the topic we refer the interested reader to the review in [2!|). What happens to large- N volume independence in 
that case ? To answer this question we find it useful to see what happens in a well defined and systematic calculation, 
and the choice of this paper is the 't Hooft model. 13 

Thus our goal in this section is to show how the volume dependent 't Hooft classical Hamiltonian Tt, derived in 
Eq. (|108[) . behaves in different cases. We first study the case of a translational invariant ground state (for any value of 
B). We do so for both the original Hamiltonian and also ask what happens if one forces the Polyakov loops to acquire 
expectation values of different sorts. We then move to discuss what happens when we allow translation symmetry to 
break. 



A. Full translation symmetry 

Since the staggered fermion Hamiltonian is invariant to translations by two lattice sites, it will be easier to discuss 
the original 'naive fermions' case in these sections. These fermions are two-component spinors and it is straight forward 
to repeat the analysis in Section IVII Bl for them. The result is that p® y becomes 2x2 matrix and that the classical 
Hamiltonian in Eq. (|119p now has the following form (here and below the trace refers to this extra 2-dimensional 



We note that the following confusion may arise : a tool one can use to make measurements in small volume theories is the Gross-Kitazawa 
'momentum feeding' trick presented in |Tol) . For example, this was used successfully in [T3| . This trick allows one to extract, from a 
zero volume theory, the meson propagator G(x) for any value of the separation x, and to measure the meson mass from the exponential 
decay in \x\. This reflects how large- N projections repackage (but not lose) the large-volume degrees of freedom into the color indices. 
A natural question now appears : can one also 'repackage', in a similar way, a baryon crystal that breaks translations into the color 
degrees of freedom of a zero volume theory ? As we shall see below the answer to this is no. This confusion arises because G(x), that 
depends on x, can be calculated from zero-volume. This, however, is a direct result of translation invariance. When the latter is broken, 
the meson propagator depends on two space-time coordinates, and the Gross-Kitazawa trick cannot be used. 
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Dirac space). 

H F /N = -gX] tr a 3P° x , x+ i + c.c. +m^2tra 1 p° x x 

X X 

+ "-T E trp«^|o:- y + Qi s |e-^-^»l. (132) 



Note that to get Eq. (|132[) we dropped the constant term that appears when one uses Eq. (|110p to write TCf in terms 
of only p (and not p) . 

Restricting to translation invariant states mean that the coherent states have 

r;, v = c% '. (133) 

Using Eqs. (|109p and (I107() this means that p® is a function of the combination x — y + Q-L s : 

p« =p(x-y + QL s ). (134) 
Combining the sums over x and Q into a single sum over the integers, we get 

-T£ tr '( r )rtOMe- e|r| . ( 135 ) 



H F /{NL S ) = tr 



Since p(r) is defined for all integer values r, we can write 

P(r) = J %e lw 9{v). 

In terms of p(p) the constraint Eq. (|B12|) . applied to naive fermions, becomes 

dp 



(136) 



trp(p) =n B + l, (137) 

Z7T 



with jiB equal to the baryon density B/L s , and the classical Hamiltonian has the form 

r27T I „2 / 



^/(7VL S ) = JJ ^tv\p{p){-a 5 Bm{j>) + ma x )]- 9 -^-j-j-^ 



g 2 N f [dp dq tr {p{p) p{q)) 



(138) 



2?r 4 sin 2 {(p-q)/2)' 

Here by^ we mean the principle value which now has a precise meaning in the form of the primed sum in Eq. (|124[) . 

The crucial point that we want to make is that Eqs. (|137p - (|138[) are independent of L s and this is the way volume 
reduction works in our model. 14 



It is easy to repeat the above steps for a gluonic state that gives nonzero expectation value to some windings 
of the Polyakov loop. This can be realized by adding a potential for the Polyakov loops, in similar lines to the 
potential suggested in [lj| (although, of course, with an opposite sign since the potential of [14| was devised to null 
all expectation values of all Polyakov loops). A potential can be chosen such that it induces, in the gluon sector, a 
spontaneous breaking of the Zjv symmetry of the form 

Z N — > 0. (139) 



We note in passing that since all physical information on the system is encoded in the classical Hamiltonian, then the theory's excitation 
spectrum will also be independent of the volume. For example, the meson spectrum, which is encoded in the l/N fluctuations around 
the minimum of 7i, will also be independent of the volume. For further details on how to extract the spectrum of mesons and glueballs 
from "H, we refer to Ref. |2- r J |. 
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In this case we have 



(P g >o = l, Vg, (140) 

and it can be easily shown that the only change this causes to Eq. (|132[) is to replace p® by the Q-independent 
function p xy , that is given by 

Qez 



Next, for a translational invariant state we write 

L 

L 



1 Ls 

P*y = -rY.P(l)^ l/L " {x - y) - (!42) 



Since p is now independent of Q, one can use Eq. (|113[l to perform the sum over Q in Eq. (|132p and one gets 

5 2 iV 1 tr (p(Z)p(fe)) 



H F /(NL S ) = T Y, tr P(0(-°"3 sin f — J + mcj-i) - 



4 L " ^ 4sin 2 (^G-fc)/L s )' 



(143) 



The important point about Eq. (| 143|) is that it depends on L s is a very strong way. Indeed, this is the reason 
Ref. [22j, which sets all Polyakov loops to be nonzero, saw a strong volume dependence of physical observables which 
was realized in a phase transition that occurs as a function of L. It is also easy to generalize this result to a breaking 
of Zm — > Zk ■ In that case L s in the l.h.s. of Eq. (| 143[) is replaced by KL S , which again means that Hp/L s depends 
strongly on L s . 

An important remark noted also in Ref. [l5j here is that when we resolved the quantum Gauss law we assumed 
that the following conditions holds 

ip a -tp b ^27rn ; n€Z (144) 

This is certainly correct for the correct solution with vanishing Polyakov loops where the field configurations that do 
not obey Eq. <| 144|) have zero Jacobian and thus zero measure (see Eq. ([76])). For nonzero Polyakov loops, however, 
Eq. (|144p can indeed be violated. This means that more zero modes (except for E J ) will be present. Since this is 
not the main topic of this paper we do not treat these here. 

To conclude, if we assume translation symmetry and unbroken Zjy, then two systems with the same ub, yet different 
volumes, will be large- N equivalent. Ref. 19] showed, in the continuum regularization, that in the translation invariant 
sector, there is a density n c B above which chiral symmetry is restored. Let us denote the dimensionless combination 
n c B / \] g 2 N by x c . 15 This can be visualized in a simple phase diagram in the space B — La/ g 2 N that we present 
in Fig. (Q]). Physics along lines of fixed slope is L independent. Nonetheless along horizontal lines, that have a 
fixed baryon number, one encounters strong L dependence realized as a phase transition. The only case where these 
two lines coincide is, of course, the B — case. We note in passing that the authors of [lj| have shown that the 
translational invariant ansatz is inconsistent and suggest that the true ground state must break translations. This 
fits well with the calculation of Ref. [18| for the B = 1 case, and we find preliminary results for B > 1 which are 
consistent with this [l6| ■ 

B. Broken translation symmetry and a crystal of baryons 

In this section we wish to emphasize that the volume independence obtained above crucially depends on the 
assumption that 

Pxy = Px-y+QL„ (145) 



Ref. [19| get x c ~ 0.0149. 
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A 

B 




L^/g^N ~ 67.3 _ 

FIG. 1: Phase diagram for chiral symmetry restoration if one assumes translation invariance. Along fixed lines of fixed ns 
there is exact volume independence. See [ly] for more details on the numeric value of the slope of the phase transition line. 

i.e. that the ground state is translation invariant. This condition is anticipated in advance from the point of view 
of orbifold projections [8( : the symmetry by which one projects QCD at large volume to QCD at small volume is 
translation invariance. 

It is useful to make a concrete example to demonstrate this point, and we first consider the case of single baryon 
in a box. At large- TV the baryon has mass of 0{N) and is a static particle 16 . As a result, its presence in the system 
spontaneously breaks translation invariance. Indeed the calculation in Ref. fl8l ] showed that the baryon wave function 
in the 't Hooft model pictorially looks like Fig. [21 where we present a sketch of the baryon density. In Ref. [l6[ we 
plan to revisit this calculation and to extend it to a finite baryon density. Let us assume that for sufficiently large 
B a crystal of baryons will form - of the sort sketched in Fig. [3] In both the single baryon and the baryon crystal 
cases the ground state wave function has a characteristic length scale A - it is the baryon width for B = 1 and the 
baryon-crystal wave-length for B > 1. These length scales must depend on g 2 N, m, and L, and it is clear that 
decreasing the box size L will change the ground state wave functions - the baryons will get squashed. This is a result 
of the compactness of space and the periodic boundary conditions. The only case where the ground state represented 
by Figs. ((2J3J) is invariant under a change in L is if A = 00 which corresponds to intact translational symmetry. 

Nonetheless, there is an interesting subtlety to the statement above : in contrast to the single baryon case, the 
baryon crystal sketched in Fig. [3] has an unbroken symmetry : shifts by A leave the ground state invariant. Thus the 
following softer form of large- N volume-independence survives in this case. 



Two systems I and II defined in volumes Li and Ln, and accommodating baryon crystals with baryon numbers Bi 
and Bji, are large-N equivalent if their baryon densities njjj — Biji/Ljjj are equal. 

The meaning of this statement is depicted pictorially in Fig. [3] : the three baryon system forms a crystal in a box of 
unit size, and so in these units A — 1/3. The statement above means we can reproduce the physics of this three-baryon 
system from a single-baryon system whose box size is 1/3 - the system bounded between the two dash-dot vertical 



We do not consider the chiral limit here, where two dimensional baryons become massless. 
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A cartoon of the baryon density for a single baryon 




x 

FIG. 2: A cartoon of the baryon density of a single baryon at large- TV where a complete breakdown of translation invariance 
takes place and large- TV volume reduction does not work at all. 



A cartoon of the baryon density for three baryons 




x 



FIG. 3: A cartoon of the baryon density of three baryons at large- TV - a partial breakdown of translation invariance : translations 
by a third of a unit still leave this state invariant. Thus large- TV reduction holds only between this system and a system with 
a single baryon on a box of size A. 

lines of Fig. [3] In these units the baryon number density of both systems is 3, and this is the only relevant parameter 
at large- N. Let us now show how this happens. 

The starting point is to modify the ansatz in Eq. (|145[) for the way p® y depends on Q, x, and y. It will be first 
useful to change coordinates and define 

r = x + QL-y, (146) 

s = x + QL + y, (147) 

Prs = P%- (148) 
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Since x, y € [1, L s ] and Q £ Z, we have r £ Z. Given one fixes the value of r, the values s obtains are r + 2, r + 4. r + 
6, . . . , r + 2L S . The periodicity in A means that 

Pr,s+A = Pr,a- (149) 

One solution to this condition is (here we assume that both L s /A and L a are even) 

Pr,s = j^T, P q {p)t VPr+lqikS - (150) 

Substituting Eq. (|150[) into the first two terms of Eq. (|132p we have 

1 st - term = j £ tr [-a 3 sm(p) (p 2A (p) - p A (p))] , (151) 

2 nd - term = J ^ tr [men (p2A(p) + Pa(p))] ■ (152) 
Substituting Eq. (|150[) into the third term of Hf gives 

tr G°g(p) PkA-q(p') 



9 =i fc=i ,/ 17 



4sin 2 ((p-p')/ 2 + 7rfc / 2 ) 



(153) 



Where we define p g <o(p) = Pq>2A{p) = 0, and po(p) = P2a(p)- The sum over A; is a result of the summation over 
x = (r + s)/2 and y = (s — r)/2 that one performs upon the substitution of Eq. QlSOp into the Coulomb interaction 
term in the Hamiltonian. Finally the global Gauss Law constraint applied on p rs becomes 

^ tr (p A {p) + PsaO)) = n B + 1. (154) 

As it stands, the classical 't Hooft Hamiltonian density Tip / {N L s ) and the Gauss law constraint, depend only on the 
parameters to, g 2 N , ns, and A - the volume is irrelevant at large- N. 

Similarly to what we saw in Section IIX Ai this volume independence is a direct result of the ansatz we took in 
Eq. (|150p that reflects a partial translation symmetry, and of the fact that all nontrivial winding of the spatial 
Polyakov loops have a zero expectation value. If we were to neglect the zero modes then all Polyakov loops would 
effectively be equal and nonzero. This will cause the infinite set of density matrices pQ y to be replaced by the Q- 
independent quantity p xy (see Eq. (|141|) 'l which depends only on the two compact coordinates x and y. As a result 
the coordinate r = x — y would have been as well, and Eq. (|150|) will be replaced by 



L 

1 

Pr.s 



TzZzZ p qtie ^Hlr/L B+qs /2A)_ (155) 
s i=l q ez 2A 

Plugging this into the classical Hamiltonian will result in a function that strongly depends on L s and this is wrong. 



C. Remark on lattice simulations of volume-reduced large-TV QCD at nonzero B or p. 

Finally, let us make a remark of relevance to lattice practitioners. Consider large- N QCD defined in a very small 
box with side L that obeys L <C 1/A.qcd- Let us also assume reduction holds for zero baryon number (this assumption 
is automatically fulfilled in two dimensions and for higher dimensions is expected to hold in modified versions of QCD 
- see [HI)- Now, force the system to accommodate a single or a few baryons. This can be done by working in the 
canonical ensemble, or by increasing the chemical potential in the grand canonical ensemble. Because the volume is 
small, the baryon density will be huge, ng ^> {^QCD) d (here d is the number of spatial dimensions), and what we 
prove in the previous section tells us that at large-JV this small volume system is equivalent to a 'standard' large- 
volume system which is extremely dense. For a lattice theory the density will be one in units of the cutoff and the 
following consequence is immediate : trying to study baryons with single-site reduced models of the Eguchi-Kawai 
type drives the theory towards the saturation regime of the lattice, where the density is of 0(1/ a). This regime is 
dominated by lattice artifacts and is unphysical. 

This may well be (partly) the reason reason why Ref. [33| , which works with a single site model at nonzero chemical 
potential p, sees either a ground state that is empty of baryons (for small p) or a ground state that is full of baryons, 
with density of 0(1/ a), and that disappears in the continuum limit. As we discuss above, such behavior is expected 
from a single-site construction, and cannot be used to study the physical regime where the baryon density is not at 
the cutoff scale. 
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X. CONCLUSIONS, SOME REMARKS, AND AN OUTLOOK 



In this paper we study the way large- N QCD depends on its volume. General arguments, such as those found in 
Ref. [|[ and in its references, tell us what are the requirements that the ground state of a large- N gauge theory needs 
in order to be volume independent. Nonetheless, we find it is useful to see how this phenomenon emerges explicitly 
in a theory which is exactly soluble. This is the reason we chose to study the 't Hooft model in this paper. 

The formalism we use is the lattice Hamiltonian formalism in axial gauge. Since we are working with a finite box 
size, the gauge fields cannot be gauged away completely and, in the gluon sector, one is left with a set of N — 1 
curvilinear quantum zero modes. These describe the N — 1 eigenvalues of the spatial Polyakov loop. Together 
with their conjugate momenta these zero modes determine the leading 0(N 2 ) dynamics. Thus a systematic large- TV 
treatment of the full Hamiltonian proceeds as follows: 

1. Treat the pure gauge Hamiltonian - find the so called 'large- N master field'. 

2. Solve for the dynamics of the fcrmions. They now interact on the background of this master field. 

Thus, the back-reaction of the fermions on the gauge fields, which is subleading in 0(1 /N) compared to the gauge 
fields, is neglected and this is consistent as long as these 1/N effects are not enhanced by any massless modes, which 
we assume to be the case. 

To perform steps (1-2) above we choose to use the lattice UV regularization and so generalized the axial gauge 
fixin g of Ref. [23| in the Hamiltonian to the lattice. We solve the fermion sector with the coherent state approach of 
Ref. 1251. The reason we choose this approach is that it is manifestly gauge invariant and easy to justify at large- iV. 
We then end up with a regularized form of the 't Hooft Hamiltonian that explicitly depends on traces of the Polyakov 
loops, and describes their interaction with quark- antiquark pairs in the background of baryons. 

Our next step was to analyze the volume dependence of the emerging 't Hooft Hamiltonian. We showed that if 
translation symmetry is intact then 

1. When the Z^ symmetry, whose order parameters are the spatial Polyakov loops, is intact, then the spatial 
coordinate in the 't Hooft Hamiltonian decompactifies and volume independence emerges. 

2. In contrast, when the Zn symmetry breaks to Zk 1 then the 't Hooft Hamiltonian has a strong volume depen- 
dence. We emphasize again that this analysis ignores a set of additional zero modes that appear in this case. 
Since this is not the main topic of this paper we do not study this issue further, but the reader should be aware 
of this point. (The focus of this paper was the case with unbroken center symmetry, which is free from this 
subtlety) 

In our case, the gauge dynamics tell us that the Polyakov loops vanish, the Zn is intact, and so volume reduction 
takes place. 

A important component in the validity of volume reduction is the fact that the ground state is translation invariant. 
In our calculation we see how this condition arises explicitly. Moreover in the case that translation invariance breaks 
down by a crystal of baryons, we show that a softer form of volume independence takes place and that, at large- N, 
instead of studying a crystal of B baryons in a volume L, one can study a single baryon in a volume L/B. In both 
cases the baryon number density is the same and together with the gauge coupling and the quark mass, these are the 
only relevant parameters determining the properties of the large- N ground state - the volume is irrelevant. 

Another aspect of large- N gauge theories which is explicitly exposed in this work is the dominance of the gauge 
fields dynamics over the fermions dynamics at large- iV, and that it also happens at nonzero baryon number. This 
means that the physics of the ground state is planar and that quark loops are suppressed. This is in contrast to the 
conjecture raised in Ref. fl7j ]. where the author suggests that quark loops are important at nonzero baryon chemical 
potential, even at N = oo. This conjecture was originally proposed to resolve a subtle apparent confusing contradiction 
between standard diagrammatic large- N arguments and the phenomenology of QCD. This confusion is absent from 
our approach to the two-dimensional case, and we see that the conjecture of Ref. [TtJ does not hold there. Briefly, 
the way this contradiction gets resolved is by non-perturbative effects, and so using perturbation theory (even if it is 
planar) is quite misleading. This means that the 'contradiction' of large- N and phenomenology is only apparent also 
in four dimensions. Further discussion on this point will be given in Ref. (20| . 



To conclude, we show how volume independence emerges in a clear and simple way in the Hamiltonian approach 
to the 't Hooft model in the lattice regularization. In the presence of a baryon crystal a partial form of volume 
independence allows one to substitute the study of the crystal of wave-length A, baryon number B, and box of size L, 
with a system of a single baryon in a box of size L = A. The latter may be useful in our companion study [1 61 ] where 
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we aim to solve the 't Hooft Hamiltonian of Eq. (| 1 08|) given a baryon number B > 1. Surprisingly, this has not been 
done yet for arbitrary quark mass (for the vicinity of the chiral limit where the two-dimensional baryons are nearly 
massless, see [lH). In fact, even in the B — 1 case, which was studied in Ref. [TH, the spatial Polyakov loops were 
set to unity. As discussed above this is inconsistent with the gauge sector dynamics and will give erroneous results at 
small enough volume (and by the equivalence mention above, at large enough densities, if one increases the density by 
fixing the baryon number and decreasing the volume). In [16| we also plan to see what will be the effect of correcting 
this issue. 

Finally we hope that knowing how two-dimensional QCD behaves at nonzero baryon number will be of value for 
studies of the physical four-dimensional system. In particular it seems that the partial independence of the QCD 
ground state on the volume, that we see emerging in the 1 + 1 case, is of a general nature and from the orbifold 
projection point of view can be anticipated on general grounds. This makes the successful modification to the Eguchi- 
Kawai reduction the four-dimensional case (of the type of [lij]) appealing, since it will allow one to study a single 
baryon in a modestly-sized box (but not of zero size) and conclude on how large- N four-dimensional QCD behaves 
at moderate/high densities. 
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APPENDIX A: HERMITICITY OF £ a AND ITS COMMUTATION RELATIONS 



In this Appendix we show that the zero mode E p L obeys 



t?" 1 T?"J 

- C/ p=0 ' rj p=0 



= 0. 



To show Eq. (|A1[) we use Eq. (|5l)) and write the difference between E ' L and its hermitian conjugate 



E 



"I 

p=0 



E 



"i 

p=0 



3 J 
"k 



±j s x x ' 1 3 L 

l T E { U U • ■ ■ <L (K Adj - w 4 K=o> 0^) 



(Al) 
(A2) 



(A3) 



To get the last line we used the fact that, within the physical Hilbert space, and except for E p = with fee [1, N— 1], all 

other components of the operator E" commute with the gauge fields (see Eqs. (|5Tj) and (pj)) ). and that E" 1 commute 
with the S operators. Finally, it is easy to show that 



(A4) 



which proves that E p L a is hermitian within the physical sector of Hilbert space. Since A^ a is real this also means that 
£ a is hermitian. 
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Let us now show that Eq. (|A2|) is obeyed within this subspace. We begin by writing the l.h.s in terms of E . The 
result is : 



^x J ^y 



(A5) 



Using the commutation relations of E' (which are the same as those of E - see Eq. (HH])) it is easy to show that the 
first term on the r.h.s. is given by 



l*-tatm=i-X;E(^ r *^"^- 



(A6) 



x k 



Here the c-number matrix is the A; th generator of SU (N) in the adjoint representation. Also, because J, J G 
[1, N — 1], we can write this term as 

l St -term = ^ E E ( A7 ) 

x k 

(fi*)„ = (Sf/i ■ T^t/t^.-.t/t^)^. (A8) 
Next we proceed to the second term and first evaluate the commutation relation of EJ with (V^^ )^ : because 

3 J 



J G [1,N - 1] then {y y ).j = (jj^ ■ ■ ■ U\& ) t , and we have 



K\ (^ Adj J = (uU ■ ■ ■ ul) fl [e';,sI] +c F ( y -i-x) [e' x \ (U, ■ ■ ■ U v -x\ 



+ C F (y - 1 - x) (Ut ■ ■ ■ U x ^) lf \e' x \ (U x ) fg \ (U x+1 ■ ■ ■ U y - X ). S] 3 



(A9) 



The first term vanishes since all E fields commute with the S fields. We now need the commutation relation between 
E and (U x ) Ad ^' . We calculated them explicitly by using Eq. ((3]), and after some algebra we find an expected result : 



'E' x \(U^) fg ] = (T i U^-) fg . 
Using this we obtain the following expression for the 2 nd term of the r.h.s. of Eq. (|A5[) 

2 nd -term = -L ^ (V x ) u (SUx ■ ■ ■ U x ^T l U x U x+1 ■ ■ ■ U y ^ x ) Ej -(!<-> J). 



x<y 



Since T{ k — —Tfc we can rearrange Eq. (jAllj) to get 



2» d -tcrm = -^E [(<®Ji (^U V* • • • U v-Jkj = "f E V*- U v~^3 E 'v 



(A10) 



(All) 



(A12) 



x<y 



x<y 



Here we got the last equality by using the fact that flu = — {£1)jt (which is a result of 2$ = -2$). 

The outcome of the above paragraphs is that the commutation relations between E p L and E p £ are proportional 
to a linear combination of (fl x ) TJ - Let us now show that (fi*)y - = if i, j G [1, iV — 1]. For that we basically write 
the definition of SI : 



(0*) tf =tr A i 5{7 1 ---r7 a: _ 1 A ff C7 t ---C7 1 t 5 t T^tr \ f ■ -U^S* tfSUi-- -U^-i 



(A13) 
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where here the trace is in the fundamental representation, and we use the fact that T% f = [A 9 ,A fc ] and the 

completeness relation of the A matrices. After some algebra we get 



(A14) 



which shows that for i,j belonging to the abelian Cartan Sub-algebra, then (0) . . = 0, and that consequently 



(A15) 



Clearly this also means that [£ a ,£b] = 0. 



APPENDIX B: CALCULATION OF Hf- 



In this section we present the calculation of TLf(Cf,'P). We begin by evaluating the expectation value of the 
numerator p a b(x)pb a (y) of Eq. (|83[) in the fermionic coherent state \Cp) = U(Cf)\0), where the reference state |0) is 
defined in Section fVII Bl For that we write 



(C F \p ab (x) Pba (y)\C F ) = {i)\U^ (C F ) Pab {x)U{C F )U^ (C F ) p ba (y) U(C F )\0). 
To calculate W{Cp) p a b{x)U(CF) we write 

u\c F )Pa b {x)u(c F ) = ~ [u\c F )r x b u{CF)u\CF)i> a x u{c F ) 



(Bl) 



Sab 

' N 



^u^{c F )r x c u{CF)u\CF)r x u{c F ) 



and use the Hadamard lemma to show that 



U\C F )V X U{C f ) = e lvJy ~ x)/L ° (e- iCF ) xv i/>. 



(B2) 



(B3) 



y ez 



This then gives 

(CF\Pab(x)Pba(y)\CF > 



E 



XV / jf> _ \ WX 



-iC F \ xv I AC 



u'ez 



iip a (v-x)/L B +iip b (x-w)/L s i\b 



V w Vv M Z^i 



N 



i<Pb(v-y)/L 3 +i(p a (y—w')/L a ^t^a _ £a& ^ p i<p c (v'->, «;')/£„,/, <=t ?/)C 



N 



|0). 



(B4) 



We proceed we need to evaluate fermionic contractions of three types: 

• Terms of type I : for the a ^ b terms of Eq. <[M|) we need to evaluate (0|V4? C V^S' |0>. Since V~|0) = ^11°) = 



for z £ S and zeSwe get 



<c#t h c io> = E w^^w 



(B5) 



Here we use the notation where x = imodI s (recall that S(S) are the set of sites that are full(empty) of 
quarks) . 

• Terms of type II : for the a — b terms of Eq. (f83|) we need to evaluate (0\ipl? ip% W^i^i |0). Here we have more 
contractions and we get 



(B6) 
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• Terms of type III : for the a = b terms of Eq. (|83|) we also need to evaluate (0|V4f V 7 " V'l/V'S' |0); a ^ c. 
Here we have only one contractions that gives 

w£ 8 <v#v£io) = E w^^w ( B7 ) 

To express Eq. (|B4|) in a compact form we define the following 'matrix densities' 

P\ 



\ = E (e lCF ) Z+pLs,a; ( e -^)^+^ +< ^ , (B8) 



%u = E (e^Y^^ ( e -^) y ' Z+{p+q)L ° . (B9) 



The hermiticity and periodicity of the matrix Cf imply that these matrix densities obey 

P q xy = S xy S q ,o ~ P q xy , (BIO) 

E W = ( BU ) 



<jez 

HEX-a 



E PL = (S + L s /2)c5 g , . (B12) 



Substituting Eqs. (|B5|) - (|B7f into Eq. (|B4|) we find that (Cf|# c |Cf) can be brought to the form 

„2 L " n q ff-q p i1nl{x-y) I L 
(Cv\HACv\ - X p -'Vo(a:-l/+«i. )/L s +i l p b (x-y+q'L„)/L s 

(C F \H C \C F ) - 4^^4rf((^,)/L s + W/is )/ 2 Xe 



+ A- yy [i n) 1 - - «'+^'| x y e M ? -o 

4i Z_, Z_, 4sin 2 (7rVi s ) r^"* J ^ 

Z=l ijij' v. / <v a 

— ^— y y ^£^4 x y e <*>.«-w _ (B13) 

4L s iV f-» ^ 4sin 2 (7r?/L s ) ^ 

To proceed we first show that the terms in the last two lines of Eq. (|B13[) are subleading. To see this note that we 
still need to perform the integral of Eq. (| 1 08|) . This will be done conveniently if we write these terms as 



n 2 N h jZ}^ (I - _L) e i2-xl/L(y-x) 

4Ls \hij 4sin 2 (7r//i s ) l Hxy ' vx FxxH 



- yy PxxPyy % x f p g p_ g , - -p„_„, ) \ . (Bi4) 

Here we have defined the fc-wound Polyakov loop operator as 

a 

The averages over the Haar measure of Pk and Pk x Pk< are known explicitly for the U (N) group [3(| [M| 

(Pk) = 4,o, (B16) 
m p \ x fx , mm(jfcl,iV) ^ 

{P k P- k >} = <5 fe fc' I 4-,o + ^2 )• ( B17 ) 

Using these results we can see that both these terms are at most of O(l). For an SU(N) group these averages are 
expected to differ by a small amount since the center of mass of the eigenvalues ip c is held fixed at zero. For 



30 



example, the average of Pk can be done explicitly : one can expand the Vandermond as a polynomial of e tiPa and see 
that (Pk} su(NX_ = if |fc| > 27V. Together with the symmetry of the Vandermond measure and the methods in 
chapter 8 of [32] we get 



(Pk) = S kfi + jj6\ k \ iN , 



(B18) 



In contrast to the two subleading terms we discussed above, the first term in Eq. (|B14|) includes a double sum over 
color indices and is thus of 0(g 2 N 2 ) ~ N. The same is true for the scaling of kinetic contribution of the fermions to 
the classical Hamiltonian : (Cf\Hk\Cf)- Using Eq. (|B3|) it is easy to show that it is given by 



(C F \H K \C F ) =N^2 

X 

which simplifies even further if we use Eq. (|B18|) 



-^Pl,x+i + cc. ) + m(-l) x p q xx 



X P n 



(B19) 



APPENDIX C: RESOLVING THE CONSTRAINTS ON p q xy 



In this appendix we wish to show how the constraints in Eqs. (|B10|[B~12|) are resolved. Our starting point is to 
write the Fourier transform 



(Cl) 



This is the most general way to express the dependence of p xv on the integer q since the range of the latter is the 
whole integers. In terms of p X y(p) the constraints of Eqs. (jB10l4BT2|) become the following infinite set of equations 
that hold for any value of p £ (0, 2tt]. 



Pxyip) = P* yx (p), 

£p*x(p) = (B + L./2), 

x=l 

^Pxy(p) Pyz(p) = Pxz(p)- 



(C2) 

(C3) 
(C4) 



The equations above can be solved by setting up, for each p, an orthogonal basis on the L s x L s space that is furnished 
by the indices x and y. We denote the orthogonal wave functions on that basis by (f> x (p), with n — 1,2,..., L s . In 
terms of these wave functions the hermitian matrices (p(p)) xy can be written as 



Pxy(p) = J2K(p)r n (p)K*(p), 



(C5) 



where r n (p) are the eigenvalues of the matrix p(p) in the basis spanned by 4>(p). Using Eq. (|C2j) we see that the 
eigenvalues r n (p) are real, while using Eq. (|C4|) . we see that they obeys 



r n (p) (r n (p)-l)=0. 
Thus r„(p) is either or 1. Finally, Eq. JC3j tells us that 

Y j r n (p) = B + L s /2, ; Vp, 

n=l 

and by ordering the cigcnfunctions 4> n (p) according to their values of r n (p) we see that 

r n (p) = 



n < B + L s /2, 
n > B + L s /2. 



(C6) 



(C7) 



(C8) 
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The end product of the discussion above is that the most general way to resolve Eqs. (|B10f[BT2|) is to write 

,2tt , B+L s /2 



J Z7T n =l 



(C9) 



with the functions forming an orthogonal basis for each value of p (note that this basis is not necessarily the 

same for each p and the precise form of </>™ (p) is determined by the dynamics of the H. minimization) . 

The last step is only done for the convenience of numerically minimizing Ti with respect to p : we discrctizc the 
momentum space p S (0, 2-7r] by using the standard form 

p dis ^ iz ° ^ with = 1,2,...,M. (CIO) 

This discretization, together with Eq. (|C9[) gives Eqs. (|121H122]) and becomes equivalent to the exact system only in 
the limit of M — > oo. 
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